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ABSTRACT 

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the 
systematically lower central oxygen abundances observed in local interacting galaxies. Here, we inves- 
tigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include 
cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to 
the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced 
inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The 
central metallicity is primarily a competition between the inflow of low-metallicity gas and enrich- 
ment from star formation. An average depression in the nuclear metallicity of ^ 0.07 is found for 
gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an 
enhancement in the central metallicity that is positively correlated with the gas content. The simula- 
tions fare reasonably well when compared to the observed mass-metallicity and separation-metallicity 
relationships, but further study is warranted. 
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1. INTRODUCTION 

The nuclear metallicities of star-forming galaxies 
are characte rized by a mas s-metall icity relatioii there- 
after, MZ iLeaueux et al.l 119791 : iRubin et al.1 Il984l : 
I2004D . ^Contributing to the scat- 
relation are interacting galaxies, 
which are consistently lower in central metallicity 
than non-inter a,cting systems of eq uivalent mass, as 
first noted by IKewlev et al.l (J2OO60 . and later con- 
firmed for ultraluminous infrared galaxi es (ULIRGs) 
JRupkc et al. 2008 1, close pairs in SPSS (lEUison et al.l 
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(jEkta fc Chengaludl2010t ). Observations such as these 
are most naturally explained as the result of vigor- 
ous merger -induced inflows of gas (Barn es fc Hernquistl 



.19911 Il996fl which rearrange the initial metallicity gra- 
dient (iShieldsl [T990I: iBellev fc Rovl IT991 iZaritskv et all 
11994 iMacArthur et al.l I2004D and "dilute" the nuclear 
metallicitv (jKewley et al.l 120061 ). In this sense, the same 
gas th a t driv e s nuclear s t arbur sts (jMihos fc Hernquist 
Il994bl 119961: llono et al.l \2004 \ and triggers central 
AGN activity and b l ack hole growth (iDi Matteo et al.l 
2005t ISprineel et al.l l2005l: IHopkins fc Hernquistl I2006r 
Hopkins et al.l 120071 ) . also results in suppressed nuclear 



metallicity. 

Simulations of merger-driven inflows of gas natu- 
rally predict a flattening of the initial metallicity gra- 
dient (jPerez et al.ll2006l : iRupke et al.llMOaHPerez et all 



|2011| ). which h as now been observed in a number of co l- 
liding systems (jKewlev et al.ll2010l:lRlIpke et al.ll2010b[ ). 
However, more than hydrodynamics are at play in deter- 
mining the nuclear metallicity evolution. To further our 
understanding, we should consider ongoing star forma- 
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tion with associated chemical enrichment, feedback from 
star formation and AGN activity, and the interchange 
of material between the stellar and gaseous phases as 
stars are born and later return material to the interstel- 
lar medium. 

The four main processes responsible for the evolu- 
tion of the nuclear metallicity are gas inflows, chemi- 
cal enrichment from star formation, gas consumption, 
and galactic outflows. These effects compete with one 
another to influence the nuclear metallicity, making it 
difficult to determine their relative importance a pri- 
ori. Numerical simulations have only recently been 
used to quantify the detailed impact of these effects on 
metallicity gr adients and nuclear m etallicities of inter- 
acting pairs. iRupke et al.l (|2010aD simulated colliding 
galaxies without chemical enrichment to explore dynam- 
ically induced changes in metallicity gradients, flnding 
that a drastic flattening can oc cur, accomp anied by a 
dro p in the nuc l ear ni etallicity. iMontuori et al.l (|2010l ) 
and iPerez et al.l (|201lD performed simulations with star 
formation and chemical enrich ment and found similar 
results to IRupke et al.l (|2010aD . but noted that the dip 
in the nuclear metallicity values can be counteracted by 
chemical enrichment from star formation. These simu- 
lations have refined our knowledge of the depressed MZ 
relation. 

In this paper, we explore the evolution of the nuclear 
metallicity during mergers using numerical simulations 
which include cooling, star formation, stellar feedback, 
and black hole growth and AGN feedback. Our approach 
allows us to systematically investigate the importance 
of gas inflows, chemical enrichment from star formation, 
gas consumption as a result of star formation, and galac- 
tic outflows. In order to unambiguously determine the 
role of metal enrichment, we have developed a stochas- 
tic method to recycle stellar particles back to the in- 
terstellar medium without requiring inter-particle mass 
mixing. The stochastic gas recycling method, which is 
designed by analogy to the widely used stochastic star 
formation method, has attributes that are distinct from 
kernel weighted mass return, making it particularly use- 



Torrey et al. 



ful for our study. 

We consider a range of situations with systemati- 
cally varied parameters to enhance our understanding of 
metallicity evolution as a natural extension of the merger 
process. We find that, for systems modeled after those 
in the local Universe, the main driver behind changes to 
the nuclear metallicity is the flow of low metallicity gas 
to the nuclear regions, or metallicity dilution. However, 
we identify two previously under-appreciated effects that 
infiuence the strength of this dilution: locking of gas and 
metals in the stellar phase, and gas and metal removal 
via stellar-driven winds. We also find that, for systems 
modeled after high redshift galaxies, the main driver of 
the nuclear metallicity shifts to chemical enrichment. In 
these simulations, the nuclear metallicity increases, con- 
trary to what is observed in the local Universe, suggest- 
ing that the interacting galaxy MZ relation may evolve 
differently at high redshifts. 

We compare our simulations directly to observations 
by synthesizing a population of progenitor galaxies with 
properties consistent with observed samples and show 
that the empirical depression in the close-pair mass- 
metallicity relation can be reproduced while account- 
ing for chemical enrichment. We find good agree- 
ment between our simulations and observations for 
both the interacting galaxy mass-metallicity relation and 
separation-metallicity relation. 

2. METHODS 

We employ a library of merger simulations carried 
out using the N-body/S mooth Particle Hydrodynamics 
(SPH) code Gadget-2 CSpringel'2005^. Gadget-2 is 
based on a formulation of SPH (Springel fc Hernguisti 
120021) which conserves both energy and entropy simul- 
taneously (when appropriate). In addition to account- 
ing for gravity and hydrodynamics, our simulations 
also include a sub-grid two-phase model of star for- 
mation, supernova feedback, radiati ve cooling of gas , 
and star formation drive n winds (jKatz et al.l 119961 : 
iSpringel fc Hernguisti [200l hereafter SH03). These fea- 
tures give a working descriptio n of the multiphase nature 
of the interstellar medium (e.g. IMcKee fc Ostrikeilll977D 
without explicitly resolving the various phases. Super- 
massive black hole sink particles are included in our sim- 
ulations and have both accretion and thermal feedback 
associated with them, but their presence does not influ- 
ence any results presented in this paper. 

The simulations presented here also include a novel 
method for returning gas to the interstellar medium 
(ISM) from (prompt) supernovae and asymptotic gi- 
ant branch (AGB) winds. While the main body of 
our simulations co de has been used for numerous other 
studie s (see, e.g., iCox et all l2006al lbl: [Robertson et al.l 
l2006allbl ia). the inclusion of time-delayed gas recycling 
is a new feature which is described. 

2.1. The Basic Model: SH03 

The primary astrophysical processes modeled in our 
simulations - including cooling, star formation, feedback, 
and galactic winds - are included as in the SH03 model. 
Because these processes are discussed extensively in that 
paper we review them only briefly here. 

The instantaneous rate of star formation for an SPH 



particle is given by 



dt 



Me 



where M* is the stellar mass, Mc is the mass of cold 
gas (i.e. some fraction of the SPH particle's mass that 
is in the cold phase as determined via the sub-grid 
model prescription of SH03), and t* is a characteristic 
star formation timescale. We adopt a value t* = 4.5 
Gyrs so that our simulation s are consistent with the 
Kennicutt-Schinidt rela tion (jSchmidti 119591: iKennicuttI 
119981 : iCox et aI]|2006bD . Star formation is assumed to 
only take place in gas which has densities above a set 
threshold, which in our case is 0.5 cm~'^. 

Equation [T] is used as the basis for a Monte Carlo 
method for actually converting SPH particles to star par- 
ticles. Each SPH particle is assigned a probability of 
turning into a star particle as follows: 



P* 



1 — cxp 



At 

tsFR 



(2) 



where At is the current simulation time step, and tsFR = 
Mgas/{dM^,/dt) is the star formation timescale. A par- 
ticle is converted if a random number drawn between 
and 1 is less than Eq. [2l 

The instantaneous rate of mass launched in star for- 
mation driven winds is given by 



dMu 
dt 



dlVL 



7- 



dt 



(3) 



where dM^,/dt is the local star formation rate and 7 is the 
mass entrainment efficiency. We typically take 7 = 0.3, 
which is the mean value for lu minous infrared galaxies 
measured bv lRupke et al.l (|2005l ) , and the resulting winds 
are ejected at Uw = 242 km/sec. 

Equation 1^ is used as the basis for a Monte Carlo 
method for wind generation. The associated probability 
of launching a particle in a wind is given by 



Pu 



1 — exp 



At 



(4) 



where t„ = Mgas I {dM^ / dt) . Particles put into winds re- 
ceive a velocity kick and are not allowed to interact with 
the surrounding medium while they leave their host en- 
vironment. In principle, gas particles launched by winds 
may still contribute to the nuclear gas phase metallicity. 
However, in practice, these particles very quickly leave 
the nuclear regions, and no longer contribute to the nu- 
clear metallicity. 

2.2. Gas Recycling 

While models for star formation and feedback are 
ubiquitous in the current generation of numerical simu- 
lations, it is much less common to include the recycling 
of mass and metals which are lost from evolving stellar 
populations, a fundamental process in the enrichment 
of the interstellar and intergalactic mediums. The few 
models that do include this "gas recycling" typically 
do so in a continuous fashion, during which, at every 
time step, a small fraction of mass and metals is 
transferred from a single stellar parti cle to a set of 
nearest SPH neighbors (e.g. .Steinmetz fc Muelleii 119941 : 



Nuclear Metallicity of Interacting Galaxies 



Kawatal [200l iKawata fc Gibsonl l2003t lOkamoto et"al 




Scannapieco et al. 120051: 



iTornatore et all l2007t 



Montuori et al.ll2010l : iPerez et al.ll2011h . This approach 
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can accurately track the spatial and temporal distribu- 
tion of mass and metals at a reasonable computational 
cost, but unfortunately, because particles can exchange 
mass with many neighbors, at many different times, it 
is impossible to trace the exact origin of any particular 
mass or metal element. 

Given a desire to unambiguously track the origin 
and evolution of interstellar metals, the work presented 
here adop ts a s tochastic approach to g as recycling (e.g. 
iLia et al.|[200llMartinez-Serrano et al.ll2008D . This less- 
used approach is designed to mirror the stochastic star 
formation scheme that is typically implemented in nu- 
merical simulations. While stochastic star formation 
probabilistically converts SPH particles into a fixed num- 
ber of collisionless stellar particles based upon the star- 
formation rate, our stochastic gas recycling converts stel- 
lar particles into SPH particles based upon theoretical 
mass return rates for a simple stellar population. When 
used in conjunction, these methods can accurately track 
the temporal and spatial interchange of baryonic material 
between stellar populations and the interstellar medium, 
a necessary element to any study of the evolution of in- 
terstellar metals. 

In our stochastic approach, the probability of turning 
an individual star particle into an SPH particle at any 
given time step is given by 



Precy 



1 — exp 



At 



recy 



where At is the current simulation time step, and 



^recy — 



dMgas/dt 



(5) 



(6) 



Here, M^,{t) is the expected stellar mass at time t, and 
dMgas/dt is the mass return rate at the same time. While 
the physical mass of the stellar particle in the simula- 
tion remains constant, the expected stellar mass identi- 
fies the amount of mass that would have been returned 
to the ISM if we implemented continuous gas recycling. 
We should make clear that there are no hybrid particles 
in the simulation, so the expected stellar mass has no 
infiuence on the simulations dynamics. It is necessary, 
however, for Eq. |6]to use the expected stellar mass in or- 
der for the stochastic recycling routine to reproduce the 
desired Starburst99 gas recycling rate in the simulations. 
In practice, the simulation selects a random num- 
ber between and 1, and performs the particle con- 
version if this number is less than Precy We calcu- 
late the mass return rate, dMgas/dt, dir ectly from Star- 
burst99 synthetic population models (iLeitherer et al.l 
[19991: IVgiauez fc Leithereii [200l ITeitherer et al.ll2010[ ). 
These models include gas recycling contributions from 
core collapse supernova and asymptotic giant branch 
(AGB) winds, but do not include contributions from 
Type la supernovae. At each time step, stellar particles 
are assigned a mass return rate determined by their age, 
metallicity, and an initial mass function (T). For the 
simulations presented here, we choose a Kroupa initial 
mass function with solar metallicity. The resulting mass 



return rate, and mass return fraction are demonstrated 
in Figure [1] 

While it will be fruitful to explore variations in the 
mass return according to varying IMF assumptions, 
real-time metallicity information, or even easily param- 
eterized functional forms of ma ss return that include 
Type la supernovae (se e, e.g., iJungwiert et al.l 120011 : 
iLeitner fc Kravtsovll201l[) we choose to allow variations 
in mass return through a simple scale factor to the return 
rate shown in Figure [TJ Specifically, we define a scalar, 
^, and employ a mass return rate normalized as such 



dM, 



gas 



c 



dt 



J7/0 fsb99{t,Z,T)dt 



/,M9(t,Z,T), (7) 



M. 



where fsb9g{t, Z,T) is the mass return rate from Star- 
burst99 as shown in Figure [T] The expected stellar mass 
in Eq. |6] is then calculated as 



M,(t) =Mo 



dMa. 



dt' 



-dt\ 



(8) 



where A/q is the original mass of the stellar particle. We 
note that for large values of C,, and for old stellar ages, 
M.^, (t) can become negative, yielding unphysical negative 
values for trecy and consequently a negative probability 
Pgas- This situation is non-catastrophic, however, be- 
cause, in practice, it leads to an immediate conversion of 
this stellar particle to SPH. 

The parameter ( allows us to control the timescale and 
efficiency of gas recycling. For example, the Kroupa IMF 
we employ dictates that 37% of the stellar mass will be 
returned to the ISM after 10 Gyr. Setting C = 0.37 
makes the normalization unity and the recycling there- 
fore tracks exactly what is shown in Fig.[T] Setting C = 1 
will produce complete (i.e., 100%) recycling of the stel- 
lar mass within 10 Gyr, while C = 100 effectively yields 
the instantaneous recycling approximation with all the 
stellar mass being returned within ^ 3 x 10^ Yr. Finally, 
and perhaps obviously, setting C — turns off recycling 
altogether giving us the fiexibility to study the impacts 
of varying amounts of recycling on metallicity evolution. 

When a star particle is returned to the gas phase, the 
particle type is instantly converted from stellar to gas at 
the end of the current time step. The fields which are 
defined for both SPH particles and star particles, such 
as the position, velocity, mass, metallicity, etc., remain 
unchanged during the conversion. For concreteness, the 
metallicity of a particle remains unchanged during the 
transition to or from the stellar state. The metallicity 
increases while in the gas phase (according to Eq|9l), but 
remains fixed while in - or making transitions to and 
from - the stellar state. 

SPH quantities must be initialized for the newly con- 
verted gas particle. In principle, we could set these fields 
to properly reflect the physical state of the gas being 
returned and explore the feedback implications that nat- 
urally result from this gas recycling model. However, in 
practice, we are resolution-limited, we already include 
feedback according to the SH03 model, and our primary 
concern is representing the overall mass-budget faith- 
fully. To this end, we give a newly formed SPH particle 
properties that will allow it to quickly homogenize into 
the surrounding medium. In particular, the entropy is 
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Figure 1. The mass recycling rate (left) and integrated mass recycling fraction (right) from a Starburst99 model for a 10 Mq stellar 
population with a Kroupa IMF. 



set equal to that of gas at a temperature of T = 50, 000 
K and a density (chosen to be the star formation crit- 
ical density) of p = 0.5 cm~^. It is important to note 
that our simulations are nearly invariant to our choice of 
these values. Changing the initial entropy by an order 
of magnitude in cither direction yields no obvious differ- 
ences in the simulations. In effect, this implementation 
of gas return serves to provide a passive source of gas for 
our evolving system, without providing strong feedback. 
During both star formation and gas recycling events, 
no particle splitting occurs. As such, the total number of 
baryon particles (stars and gas) is conserved throughout 
the simulation. Effectively, baryon particles are permit- 
ted to flip back and forth between the gas and stellar 
phases, according to the star formation and gas recycling 
rates. Because baryon particle number is conserved and 
no mass is transferred between particles, mapping a par- 
ticle's initial position to a final position a trivial task, 
regardless of the number of times it was it was turned 
into a star particle or returned to the gas phase. Using 
this stochastic gas recycling method provides the distinct 
advantage that we can retrace all mass to a unique initial 
position regardless of the star formation or gas recycling 
history. 

2.3. Metallicity Enrichment 

We implement a method for calculating the metal en- 
richment based on the instantaneous star formation rate 
to determine the metallicity enrichment rate. Unlike star 
formation, galactic winds, and gas return, metal enrich- 
ment is carried out in a continuous fashion, where the 
metal formation rate is given by 



dMz 
dt 



dM, 



y- 



dt 



(9) 



and y is the metal yield and dM,,/ dt is the instantaneous 
star formation rate. For our simulations, we use a fixed 
yield of j/ = 0.02. The metallicity will increase wherever 
there is ongoing star formation, as determined via the 



SII03 star formation model. As such, diffuse gas will not 
be star forming, and will therefore not enrich (see SII03 
for details). The scalar metallicity is updated at each 
time step using an Eulerian integration scheme. This 
scalar metallicity does not track independent species, but 
instead provides a single metallicity value proportional 
to the integrated star formation rate. This metal enrich- 
ment scheme is strictly independent of our gas recycling 
routine. 

The metallicity enrichment does not affect the evolu- 
tion of the simulation because there are no metallicity 
dependent dynamical processes included, such as metal 
line cooling. Because no mixing is allowed, nearby par- 
ticles may have large metallicity variations. This is not 
problematic, as long as we interpret the metallicity at a 
given location to be a kernel weighted average of nearby 
particles, as is traditional for determining fluid quantities 
in SPH. Mixing would homogenize the individual parti- 
cle metallicities, while leaving the average value of the 
kernel weighted metallicity unchanged. 

Using this enrichment scheme, we are able to cleanly 
decompose the metallicity of an SPH particle into con- 
tributions from its initial metallicity and star formation 
history. More important, we are able to scale our metal 
yield and modify our initial metallicity setup in our post- 
processing analysis, without requiring additional simula- 
tions. This is possible because all of the gas content of 
an SPH particle has an unambiguous and unique initial 
location and star formation history, which is not true 
when particle mass mixing is used. In later sections, 
our ability to arbitrarily modify initial metallicity gradi- 
ents of our progenitor galaxies is critical in allowing us 
to thoroughly sample the range of metallicity gradients 
of our progenitor galaxies, without large computational 
requirements. 

3. ISOLATED GALAXIES 

In our analysis, we take the nuclear region of a galaxy 
to be a sphere of radius 1 kpc about the galactic center. 
We determine the nuclear metallicity from the star for- 
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mation rate weighted average of all gas particles inside 
this sphere. A 1 kpc spherical r egion is used f or con - 
sistency with the observations of iKewlev et alj ()2006[ ) . 
However, it should also be noted that our results would 
not fundamentally change for slightly larger or smaller 
definitions of the nuclear region. 

We use star formation rate weighted averages to 
mimic observations of HII regions, which naturally se- 
lect star forming gas. Unless otherwise stated, all nu- 
clear metallicities quoted in this paper are gas-phase 
and star formation rate weighted. In particular, the 
central depressions in gas-phase metallicity presented 
here should not be confused with enhancements in stel- 
lar metallicity seen in late-stage mergers and relaxed 
elliptical galaxies. These increases in stellar metallic- 
ity are relics from the merger-driven starbursts that 
leave behind centra l stellar cu sps in merger r emnants 
(|Mihos fc Hernauisti 1994a: Hopkms et al.112001 and el- 
hptical galaxies (jHopkins et al.l l2009al lcf) . These star- 
bursts occur at late stages in a merger, following coa- 
lescence, from enriched, star-forming gas, yielding a cen- 
tral population of young, metal-enhanced stars (see, e.g. 
Figs. 27 and 28 in Hopkins et al. 2009a). 

The star formation rate weighted nuclear metallicity is 
defined as 

^_ J^p.ir)Zir)dV ^^^^ 

Iv P* (^ dV 

where p^{x,y,z) is the star formation rate density, 
Z (x, y, z) is the fraction of gas-mass in metals, and the 
integral is performed over the 1 kpc spherical region 
about the galaxy's center. The star formation rate den- 
sity is calculated as 



p*{r) 



^—' dt 



\r-rj\,hj) 



(U) 



where W {\f — rj\ , hj) is a smoothing kernel function, h 
is a smoothing length, and dM^,/dt is a particle's instan- 
taneous star formation rate (discussed in §2.ip . For our 
analysis, we use a galaxy's central black hole to define the 
galactic center; however, our results are unchanged if we 
use the potential minimum, or any other reasonable mea- 
sure of the galactic center. This approach neglects pro- 
jection effects, but instead provides information about 
the "true" nuclear region. 

The metallicity of each SPH particle consists of two 
separate contributions: an initial metallicity and an en- 
riched metallicity. The initial metallicity of an SPH par- 
ticle is completely determined by the particle's initial 
radial position (discussed in §3.ip and is unchanging in 
time, while the enriched metallicity grows owing to star 
formation. The metallicity of an SPH particle as a func- 
tion of time is given by 



Zgas (t) — Zinit + 



y 



M, 



gas Jo 



•^.. (., 



where d]\'U{t')/dt' is the particle's complete time depen- 
dent star formation rate history. The integration of the 
rightmost term in Eq. [T^ should be carried out start- 
ing at the beginning of the simulation (i.e. t — 0), re- 
gardless of whether a particle has been recycled or not. 
Since dM^(t')/dt' is non-zero only when a particle is ac- 



tively star forming (i.e. in the gas state and above the 
star formation threshold density), periods of time when a 
particle is in the stellar state or below the star formation 
threshold density (see ij2.ip will not contribute to Eq. \T% 
The metallicity of the nuclear region is calculated sim- 
ilarly to the metallicity of an individual SPH particle, 
with an additional sum over an ensemble of particles 
weighted by their star formation rates. We explicitly 
break the nuclear metallicity into two components 



Zj — ^^ \ -t- 



E.M* 



E.M,,, 



(13) 



where the sum is performed over the particles that fall 
within the nuclear region. The first term contains the 
contribution from the initial metallicity gradient of our 
progenitor galaxies while the second term contains the 
contribution from star formation induced metallicity en- 
richment. For clarity throughout, we call the first term 
the "dynamical metallicity" , as it is a dynamical result 
of the initially assumed metallicity gradient, while we 
call the second term the "enriched metallicity", as it is 
a product of metal enrichment. For concreteness, the 
metallicity gradients for one of our "C" isolated galaxies 
evolved in isolation for 2 Gyrs is shown in Figure [5] with 
contributions from the dynamical and enriched compo- 
nents explicitly shown (see ij3.2p . 

3.1. Isolated Galaxy Setup 

The isolated galaxies used ii i this paper are modeled 
following the analytical work oflMo et al.l (Il998f) em plov- 
ing the procedure outlined in lSpringel et al.l (|2005[ ). Our 
fiducial galaxy consists of a dark matter halo, an em- 
bedded rotationally supported exponential disk, a stellar 
bulge, and a central supermassive black hole. 

We construct a set of four isolated galaxies with total 
system masses ranging from '^ 10^^ Mq to ~ 1O^^M0, as 
outlined in Table [TJ Although varied in mass, all sys- 
tems are constructed to be self-similar in order to isolate 
the effects of mass from other quantities that may cor- 
relate with mass in observed systems. The total mass 
of the disk (stars and gas combined) is chosen to be a 
constant fraction (4.0%) of the halo mass for all progeni- 
tor galaxies, with four settings for the initial gas fraction 
(8%, 20%, 40%, and 60%). All halos have spin param- 
eters equal to A = 0.05, which effectively sets the disk 
radial scale length. The initial stellar disk scale height 
is set to a fixed fraction (0.2) of the initial disk scale 
length, while the gaseous scale he ight is determin ed by 
satisfying hydrostatic equilibrium (jSpringel et al.l HOOS). 
Stellar bulges are included as part of the fiducial galaxy 
setup and are given a fixed fraction (1%) of the system 
mass. The resulting disk setup has a Toomre Q parame- 
ter that varies as a function of galactocentric radius, but 
is everywhere greater than unity - indicating stability 
against axisymmetric perturba tions. 

The dark matter halos follow iHernguistJ ()1990f ) profiles 
with concentration indices of c = 10. The stellar and 
gaseous disk components are modeled with exponential 
surface density profiles 



(r) 



M, 



2TTh 



2^^P[-h 



(14) 



where M*^g is the total mass of stars and gas in the 
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Table 1 

Progenitor Disk Properties 



Disk Identifier 



A 
B 
C 
D 

C2 
C3 

C4 



Total Halo 

Mass [Mq] 



Initial Stellar 

Mass [Mq] 



Initial Disk 
Gas Fraction 



h 
[kpc] 



Nr, 



Na 



N, 



2.3 X 10" 
5.1 X 10" 
9.5 X 10" 
13.5 X 10" 

9.5 X 10" 
9.5 X 10" 
9.5 X 10" 



1.1 X 10"' 
2.4 X 10^° 
4.4 X 10^° 

6.2 X 10^0 

4.0 X 10"' 

3.2 X 10^0 

2.4 X 10^0 



10% 
10% 
10% 
10% 

20% 
40% 
60% 



2.4 532,500 30,000 145,500 

3.2 532,500 30,000 145,500 

3.9 532,500 30,000 145,500 

4.4 532,500 30,000 145,500 

3.9 532,500 40,000 135,500 

3.9 532,500 60,000 115,500 

3.9 532,500 80,000 95,500 



disk, respectively, and h is their com mon scale length 
set by the disk's angul ar momentum (jMo et al.l 119981 : 
iSpringel fc W hite" 1999^. Stellar bulges are taken to be 
spherical Hcrncniist (1990) profiles, where the bulge 
scale length is given in terms of the halo scale length, 
Qb = 0.2a, and the bulge mass fraction to be 1% of the 
halo mass. 

The initial metallicity profile of our progenitor galaxies 
enters as an assumption. We model the metallicity profile 
of the disk as an exponential 



(15) 



Z{r) = Zq exp ( - — 



where the central metallicity, Zq, and the metallicity 
scale length, h^, are chosen to be consistent with obser- 
vations. Isolated galaxies show dispersion in their radial 
metallicity gradients. We could include this dispersion 
in our simulations by adding in a Gaussian error term to 
Eq. 1151 However, no dispersion is used when setting the 
initial metallicity profile for our isolated galaxies because 
this would add an extra element of uncertainty into our 
models and add noise to our results, without yielding 
systematic changes or further insight. 

We assume that Oxygen makes up 30% of the metals by 
mass. The constant Oxygen mass fraction is a necessary 
assumption of our models, because we track only one 
scalar metallicity value. We note that our results are 
not particularly sensitive to our choice of the Oxygen 
mass fraction, since this is merely a normalization of the 
overall metallicity. 

We fix the central metallicity values to the observed re- 
lation of lTremonti et all ()2004D so that the central metal- 
licity is given in terms of the stellar mass 

12 -I- log (O/H) = -1.492 + 1.847X - 0.08026x2 (16) 

where x = logiQM^,. Initially, we use the exact value 
given by equation ([TS]) when discussing the process of 
nuclear metallicity evolution in ^ However, in [J5] we 
pick the central metallicity of each progenitor galaxy 
from a Gaussian distribution with a standard devia- 
tion, (7 = 0.1 dex, taken from iTremonti et al.l (I2004D. 
Simila rly, we use the average result from iZaritskv et all 
(|1994D of h, = h/0.2 in g] However, in glwe pick the 
slope of the metallicity gradient from a Gaussian distri- 
bution with a standard deviation of ct = O.S/ij. 



3.2. Isolated Galaxy Evolution 

While the metallicity evolution of isolated galaxies is 
not the focal point of this paper, we take a moment in 
this section to demonstrate that our galaxy models are 
indeed stable, and evolve very little in isolation. We 
specifically demonstrate that the galaxy models do not 
develop strong bars since these structures can lead to 
significant radial infiows of gas. Such processes obscure 
the true origin of any central metallicity evolution, and 
what is specifically a result of the merger. 

The evolved gas surface density, seen in Figure |3l 
demonstrates that our disks are both stable and free of 
strong bars over at least 2 Gyrs, the typical duration 
of the galaxy major mergers we study. As a result, the 
metallicity gradients of our isolated systems evolve very 
little, as seen in Figure [2j We further demonstrate the 
characteristics of our isolated disks in Figure |4] where we 
show the star formation rates and gas fractions for the 
C2, C3, and C4 systems. 

To understand where the small amount of metallicity 
evolution does come from. Figure [5] also shows the evolu- 
tion of the initial metallicity gradient, or the "dynamical" 
metallicity (middle panel), and the contribution from on- 
going star formation, or the "enriched" metallicity (right 
panel). From this information it is easy to see that the 
small increase in central metallicity and gradient is a re- 
sult of chemical enrichment from ongoing star formation, 
while the initial metallicity gradient remains nearly un- 
changed. 

Motivation to emphasize the stability of our initial disk 
models stems from the fact that the dis ks used in pre- 
vious studies (e.g. iMontuori et al.l l2010f l appear to be 
bar unstable over short timescales (see their Figure 2 in 
their Appendix). These bars result in strong radial mix- 
ing and a flattening of their radial metallicity profiles 
(see Figure 1 in their Appendix). While bar induced gas 
infiows are indeed a valid physical mechanisms for mod- 
ifying the nuclear metallicity, we wish to eliminate this 
complication from our study to make our results easier 
to interpret. 

The evolution of the nuclear metallicity of the four pro- 
genitor disks, described in Table [H are demonstrated in 
Figure [5l The nuclear metallicity of all isolated disks is 
well described by a mild and monotonic increase over at 
least 2 Gyrs. The steady behavior of the isolated galaxy 
nuclear metallicity ensures that any strong changes dur- 



Nuclear Metallicity of Interacting Galaxies 




4 6 8 10 
R (kpc) 



4 6 8 
R (kpc) 



4 6 8 10 12 
R (kpc) 



Figure 2. Metallicity profiles for the "C" isolated galaxies evolved for 2 Gyr. The total (left), dynamical (center), and enriched (right) 
metallicities are shown separately. The horizontal dashed line denotes the mean initial metallicity for the disk. The stability of the 
dynamical metallicity indicates that radial gas flows are minimal in our isolated systems, while the total metallicity gradient increases 
over time owing to chemical enrichment. A clear distinction is made between the dynamical metallicity, which is a product of our initial 
conditions, from the enriched metallicity, which is a product of ongoing star formation in the simulation. 



ing the merger simulations are a product of the merger 
process, rather than disk instability. 

While exploring the behavior of isolated disks, we can 
also demonstrate the effects and capabilities of stochastic 
gas recycling. When evolving the "C" systems (i.e. C2, 
C3, and C4 from Table [T]) we vary the gas recycling pa- 
rameter and examine the resulting gas fraction and star 
formation rate evolution. The results, shown in Figure|4l 
demonstrate that the gas is locked into stars more effi- 
ciently when the gas recycling parameter is small. As a 
result, the star formation rates, which largely depend on 
the amount of gas available for star formation, are con- 
sistently larger for increasing values of the gas recycling 
parameter. Simulations that use efficient gas recycling 
maintain larger gas fractions and star formation rates 
throughout the simulation. 

In the limit of very large gas recycling parameters, all 
newly formed stellar material is quickly returned to the 
ISM. As a result, the gas fraction for large gas recycling 
parameters stays nearly constant. The corresponding 
star formation rates slowly increase as gas naturally cools 
and condenses over time. 

4. MERGING GALAXIES 

While keeping all of the conventions for nuclear metal- 
licity defined in Sj3l we now consider the evolution of 
galactic nuclear metallicity for merging systems. Since 
our interacting systems are spatially extended and over- 
lap substantially during the interaction, they do not fol- 
low Keplerian trajectories. However, we use the termi- 
nology of Keplerian orbits to clearly describe the merger 
setup. 

To completely specify our merging setup, we must 
specify the properties of the merger orbit (2 parameters), 
and the relative orientation of each galaxy with respect 
to the orbit (2 parameter for each galaxy). We start by 
specifying that the systems will be on zero energy orbits 
(i.e. eccentricity value of unity for Keplerian objects). 
The orbital angular momentum is then set by picking 
a value for the impact parameter, assuming Keplerian 
trajectories (the "real" impact parameter found in the 
simulations will be larger than this value). We hold the 
impact parameter fixed (5 kpc) for all simulations. The 
spin angular momentum vector of each galaxy is varied 



with respect to the orbital angular momentum vector ac- 
cording to the orientations detailed in Table [2] as visually 
depicted in Figure [51 

The parameter space of mergers is quite large. To ex- 
plore this parameter space, without using an excessively 
large number of simulations, we first perform a detailed 
exploration of a single case. We then systematically vary 
merger parameters, such as the disk orientation, progen- 
itor disk mass, merger mass ratio, orbital angular mo- 
mentum, and so forth. While this does not necessarily 
cover all of merger parameter space, it does enable us to 
paint a coherent picture detailing the relation between 
the evolution of the merger and the evolution of the nu- 
clear metallicity. 

We select our fiducial galaxy merger to be two identical 
"B" disks merging on the "e" trajectory (see Tables [1] 
and[2|). This setup is neither average nor special, and is 
simply selected to demonstrate in detail the relationship 
between the merger state and evolution of the nuclear 
metallicity. While HA.1\ is restricted to understanding a 
single merger evolution, i j4.2ll4.4l show that much of the 
nuclear metallicity evolution can be understood in terms 
of generic merger properties that will continue to drive 
the nuclear metallicity evolution as we move beyond this 
fiducial setup. 

4.1. Metallicity Evolution in Merging Systems 

In stark contrast to the isolated systems presented in 
i j3.21 the nuclear metallicity of interacting galaxies is 
a complicated function of time, as shown in Figure [71 
Moreover, there is a clear change in the evolutionary be- 
havior of the nuclear metallicity following pericenter pas- 
sage when the disks stop behaving like isolated systems, 
and enter a period of evolution dominated by the merger 
dynamics. 

Before delving further into the merger details, we note 
that there is a characteristic "double dip" shape to the 
nuclear metallicity evolution. This is caused by the ongo- 
ing competition between metallicity dilution and chem- 
ical enrichment. While the exact evolutionary track 
shown here is specific to this particular merger setup, 
the characteristic double dip in the metallicity will ap- 
pear repeatedly in many other cases. In general, galactic 
close passages lead to gas inflows (yielding inflows of low 
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Figure 3. The gas surface density is sliown for our isolated "C" disk evolved for 2 Gyrs. The disks evolve stably and without bar formation 
or significant radial gas inflow for 2 Gyrs. 
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Figure 4. Star formation rates and gas fractions for the C2, C3, and C4 systems evolved for 2 Gyrs with varying gas recycling parameters. 
The gas return parameters, identified in the legend, are varied from (^ = 0.0 (i.e. no gas recycling) to f = 100.0 (i.e. nearly instantaneous 
recycling). Unless otherwise noted, subsequent sections and simulations use (^ = 0.3, which is the unsealed value taken from the Starburst99 
simulations. 



metallicity gas) which are quickly followed by star for- 
mation (yielding chemical enrichment). 

During the period of time between first pericenter pas- 
sage and final coalescence, the nuclear metallicity evolu- 
tion is driven by: 1) Radial inflows of low metallicity gas 
caused by the tidal interaction, 2) Chemical enrichment 
caused by star formation activity, 3) Locking of metals 
and gas in the stellar phase, and 4) Removal of gas and 
metals via galactic winds. The combined effect of these 
four mechanisms, which are discussed in detail in the fol- 
lowing sub-sections, determines the relative importance 
of metallicity dilution and chemical enrichment, which 
ultimately decide when and where the nuclear metallic- 
ity will be diluted or enhanced. 

4.1.1. Metallicity Dilution and Chemical Enrichment 

The influence of radial gas inflows and chemical en- 
richment are demonstrated in Figure |51 which shows the 
galactic separation, rate of change of nuclear metallicity, 
nuclear star formation rate, and nuclear gas inflow rate. 



In terms of basic merger properties, there is an influx 
of gas following first pericenter passage and periods of 
gas inflow associated with each additional close passage. 
These gas inflows give rise to high nuclear star formation 
rates. While these previous points have been studied 
extensively in other pape rs (Barnes & Hernquist I1991L 
[Tool iMihos fc Hernauistl[T994b, .1996; Jono et al.l|200l . 



we instead focus here on the influence that these generic 
merger properties have on the evolution of the nuclear 
metallicity. Specifically, times of strong gas inflow corre- 
spond to periods of nuclear metallicity depression, while 
high star formation activity aligns with nuclear metallic- 
ity enhancement. These qualitative relationships remain 
true as the merger parameters are varied. 

Previous studies have found that the depression in the 
nuclear metallicity is correlated w ith the ma ss of gas 
that migrated to the nuclear region (jRupke et a l. 2010a). 
This result is reproduced in our simulations when we 
neglect star formation (similar to the red line in Fig- 
ure [7]) and is a clear-cut demonstration of metallicity 
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Figure 5. Nuclear metallicities for the A, B, C, and D, isolated 
galaxies evolved for 2 Gyrs is shown. The nuclear metallicity of our 
isolated model galaxies increases monotonically with time owing to 
chemical enrichment, as demonstrated in Figure [5] 




Figure 6. Schematic representation of the merger setup. Our 
coordinate system is defined by the plane of the merger orbit, and 
each disk's orientation is independently adjusted with respect to 
this plane. 



Table 2 

Orbital Orientations 

Orientation 8i (pi 82 (^2 

Identifier [dog] [dog] [dcg] [deg] 
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dilution. However, when we also consider contributions 
from chemical enrichment and the consumption of gas 
via star formation, this correlation disappears. 

Inspection of Figure[5]shows that, in general, periods of 
strong gas inflow occur simultaneously with times of on- 
going nuclear metallicity depression. Similarly, periods 
of high star formation give rise to enhancements in the 
nuclear metallicity. Interestingly, neither the nuclear star 
formation rate nor the nuclear mass inflow rate correlate 
well with the rate of change of the nuclear metallicity 
because they exchange roles playing the dominate driver 
of the nuclear metallicity throughout the merger. The 
reason is that metallicity dilution, chemical enrichment, 
locking of gas in the stellar phase, and galactic winds 




-0.4 -0.2 

Time (Gyrs) 

Figure 7. The total (black), dynamical (red), and enriched (blue) 
nuclear metallicities as a function of time for the "BBe" merger 
simulation. The two disks, which do not have identical metallic- 
ity evolutionary tracks, are distinguished by the solid and dashed 
lines. First pericenter passage and final coalescence are marked 
with vertical dashed and dot-dashed lines, respectively. 



all influence metallicity evolution. While more involved 
parameters (e.g. the difference between the star forma- 
tion and nuclear mass inflow rate) yield better, albeit 
imperfect, correlations, these trends can be misleading 
and they tend to over-complicate a fairly simple point. 
In particular, changes to the nuclear metallicity at any 
time during the interaction can be understood by ex- 
amining the nuclear star formation and gas inflow rates. 
Furthermore, the nuclear star formation and gas inflow 
rates are naturally explained via the well studied merger 
process. 

From this single merger example, it seems plausible 
that the role of metallicity dilution and chemical enrich- 
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Figure 8. Several diagnostics for assessing the nietallicity evolution of our fiducial merging system are shown. The top row shows contour 
plots of the gas density, with lines indicating the stage of the merger. From top to bottom, the subsequent time series show the galactic 
nuclear separation, the rate of change of the nuclear metallicity, nuclear star formation rate, and nuclear gas inflow rate. Pericenter passage 
and final coalescence are denoted by dashed and dot-dashed lines, respectively. In general, periods of ongoing nuclear metallicity dilution 
can be associated with strong nuclear gas inflows, while times of ongoing nuclear metallicity enhancement are associated with high nuclear 
star formation rates. 



ment are natural consequences of the merger process. 
However, by inspection of Figure[71 neither effect is over- 
whelmingly dominant. Hence, in the following analysis 
we pay attention not only to metallicity dilution and 
chemical enrichment, but also effects such as the lock- 
ing of gas in stars and the launching of galactic winds 
that work to modulate the efficiency of these processes. 
Instead of searching for correlations between the changes 
in the metallicity and these distinct processes, we vary 



parameters in our simulations to test their influence on 
the evolution of the nuclear metallicity. 

4.1.2. Gas Consumption 

Any process that can modify the gas reservoir in the 
nuclear region can affect the metallicity measurement. 
Here, we consider the influence of locking gas in the stel- 
lar phase when determining the evolution of the nuclear 
metallicity. Star formation lowers the reservoir of gas 
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in the nuclear region, and hence amplifies the effect of 
metallicity dilution that occurs when low metallicity gas 
floods the central region. While this is definitely a phys- 
ical effect, its magnitude will be misjudged if the simu- 
lations do not properly account for gas recycling. Using 
the stochastic gas recycling method outlined in §2.2) we 
are able to modulate the efficiency with which gas be- 
comes locked in the stellar phase without changing our 
star formation efficiency or feedback prescriptions. In- 
stead, we change the gas recycling parameter, allowing 
gas to be returned from the stellar state back to the ISM. 

We first consider a simulation where newly formed stel- 
lar material is instantaneously returned to the ISM, such 
that no gas will be trapped in the stellar phase. At the 
opposite extreme, we perform a simulation with star for- 
mation, but no gas recycling (i.e. C — 0.0), which would 
trap the largest amount of gas in the stellar phase. In be- 
tween these extremes, we consider three values for the gas 
recycling fraction, ^, that show how the nuclear metal- 
licity evolution changes as the efficiency with which gas 
is locked in the stellar state is varied. 

When the gas recycling parameter is set to very large 
values, the resulting enriched metallicity becomes un- 
physically high. This occurs when many generations 
of star formation are allowed to occur over very short 
timescales as the gas is quickly recycled into the ISM 
repeatedly. However, the resulting dynamical metallic- 
ity remains physical and meaningful, as it describes the 
change in the nuclear metallicity that will occur as a 
function of the efficiency with which material is locked in 
the stellar phase. In fact, dynamical metallicity given by 
our simulation with instantaneous gas recycling (see the 
black line in th e center panel of Figu re [9]) is very similar 
to the results of lRupke et al.l (|2010al ). These simulations 
capture changes to the nuclear metallicity caused by di- 
lution without any contributions or contamination from 
the creation of stars. In this case, the nuclear metallicity 
is well-correlated with the mas s of gas in the nuclear re- 
gion, as noted bv' Rupke et al.l (|2010al see their Fig. 2). 
This indicates that the initially high nuclear metallicity 
is diluted by inflows of low metallicity gas. However, as 
the gas recycling parameter is dialed back to lower val- 
ues, gas is more efficiently trapped in the stellar phase 
which increases the effect of metallicity dilution. 

The simulation with star formation and no gas recy- 
cling brackets the upper end for the amount of gas that 
can be trapped in the stellar phase. In this case, the 
dynamical metallicity (red dotted line in the center plot 
of Figure ^ is lower than the dynamical metallicities 
of all other simulations and falls ^-^0.2 dex below that 
of the nuclear metallicity with instantaneous recycling. 
This difference is caused purely by the effect of locking 
gas in the stellar state. When no gas recycling is used, 
a substantial fraction of the inner gaseous reservoir is 
converted into the stellar state, allowing influx of low 
metallicity gas to have stronger influences. In addition, 
the enriched metallicity is lower than any other simula- 
tion because the resulting star formation rates are lower 
when gas is efficiently locked in the stellar state. The 
resulting total metallicity is below the total metallicity 
where gas recycling is used. 

As gas recycling is included at intermediate values, 
both the dynamical metallicity and enriched metallic- 
ity increase. The increase in the dynamical metallicity 



for larger gas return fractions is well-motivated by the 
simple physical arguments presented above. Specifically, 
increasing the gas recycling fraction lowers the efficiency 
with gas is depleted from the nuclear region into the 
stellar state. Hence, as the gas recycling fraction is in- 
creased, the mass of gas in the inner reservoir increases, 
and the influence of metallicity dilution is decreased. 

The increase in the enriched metallicity for larger gas 
return fractions is similarly explained. Gas that resides 
in the nuclear region will tend to remain there for longer 
periods of time (on average) and achieve higher instanta- 
neous and integrated star formation rates (hence, higher 
metallicities). While this sheds light on the influence of 
gas recycling, we would like to use physically motivated 
values for the gas recycling fraction in what follows. We 
take the unsealed Starburst99 mass return fraction as 
the fiducial mass return fraction used for the remainder 
of this paper. 

4.1.3. Galactic Winds 

Galactic winds are thought to play a leading role in 
the evolution of the distribution of metals throughout 
cosmic time. Starburst-driven, or AGN-driven winds 
are a ubiquitous phenomena that act to transport met- 
als from the depths of galaxies into the interstellar , or 
intergalactic medium (see, e.g.. iVeilleux et aLll2005l for 
a comprehensive review). Galactic winds may also play 
a prominent role in shaping the m ass-metallicity rela- 
tionship (see, e.g.. lDave et a l. 2006) and because galaxy 
mergers induce episodes of enhanced star formation, they 
may also induce periods of significant feedback and galac- 
tic outflows thereby altering the central metallicity evo- 
lution of galaxy mergers. 

In the context of this paper, galactic winds can remove 
gas from central regions and modify the rese rvoir of ma- 
terial contributing to the nuclear metallicity (jPerez et al.l 
1201 If ). In this capacity, the influence of galactic winds is 
identical to the locking of gas in the stellar state. To 
demonstrate the impact of galactic winds on the nuclear 
metallicity, we present a set of simulations where we vary 
the mass entrainment, 7. In these simulations we do in- 
clude gas recycling at the flducial level, but note that the 
our conclusions are unchanged if this effect is turned off. 

The dynamical metallicity is not substantially affected 
by galactic winds. Since star formation is already fairly 
efficient at locking gas in the stellar state, galactic winds 
provide only an incremental modification to this effect. 
However, notable changes occur in the enriched metallic- 
ity. Specifically, the enriched metallicity is depressed for 
large values of the mass entrainment because gas with 
high star formation rates is either quickly converted into 
the stellar phase or launched in a wind. Therefore, the 
integrated star formation rates and enriched metallicity 
values for gas particles remains lower, on average, for 
higher mass entrainment values. 

The resulting nuclear metallicity, shown in Figure lTUl is 
depressed as mass entrainment is increased. This change 
is primarily caused by changes to the enriched metallic- 
ity. However, there is one additional effect that distin- 
guishes galactic winds from locking of material in the 
stellar phase. While gas that is permanently locked in 
the stellar state will forever be unable to contribute to 
the nuclear metallicity or star formation, gas that has 
been launched in a wind will ultimately fall back onto 
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Figure 9. The total (left), dynamical (center), and enriched (right) nuclear metallicities as a function of time for "BBe" merger simulations 
arc shown. The various lines show the nuclear metallicity evolution with several values of the gas recycling coefficient as well as complete 
and instantaneous gas recycling. For very large values of the gas recycling coefficient (i.e. f >> 1), stellar particles are returned to the 
gas phase quickly. As the gas recycling coefficient is increased, the results approach the case of instantaneous and complete gas recycling. 
First pericenter passage and final coalescence are marked with vertical dashed and dot-dashed lines, respectively. 
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Figure 10. The nuclear metallicities are shown for varied wind launching efficiencies. First pericenter passage and final coalescence are 
marked with vertical dashed and dot-dashed lines, respectively. 



the disk, or even into the nuclear region. We see this 
effect having an influence in the enriched metallicities of 
the highest mass entrainment value simulations at late 
times. Specifically, as gas inflow from previously ejected 
wind material falls upon the nuclear region, the star for- 
mation rates and enriched metallicities are pushed to- 
ward higher values. 

It should be noted that the galactic wind launching 
prescription used in this work does not describe a situ- 
ation where the galactic wind predominately ejects en- 
riched gas. In other words, the metallicity of the wind is 
the same as the ambient metallicity of the star forming 
region from which the wind was launched. This pre- 
scription could underestimate the efficiency with which 
metals are pushed out of a galaxy and into the intergalac- 
tic medium. However, this wind launching prescription 
is straightforward to understand and has the direct and 
notable effect of removing gas from the nuclear region. 

Galactic winds do little to change the radial rearrange- 
ment of gas which attends a galaxy merger, as shown by 
the relative insensitivity of the dynamical metallicity to 
the efficiency of the winds. The only appreciable effects 



are on the enriched metallicity. In the following anal- 
ysis, we take 7 = 0.3 as the fiducial mass entrainment 
value (|Rupke et al.ll2005D . 

4.2. Disk Spin Orientation 

While the previous analysis focused on a single merger 
setup, the effect of merger orientation can be tested us- 
ing a set of identical "B" galaxies merging from 16 ori- 
entations (a-p, as detailed in Table [2]). All orientations 
maintain unique metallicity evolutions, as the specifics 
governing the magnitude and timing of the gas inflows 
and starbursts are orientation dependent. Here we will 
consider which aspects of our fiducial merger's metallic- 
ity evolution are found to hold for various merger ori- 
entations. To test this, we first examine the metallicity 
evolution for the four orientations shown in Figure 1111 

As is evident from this figure, three of the four disks 
show strong depressions in their nuclear metallicity fol- 
lowing pericentric passage. This is caused by the strong 
tidal interaction and resulting gas inflows which are a 
natural byproduct of the merger process. The disk not 
experiencing a post pericentric nuclear metallicity de- 
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Figure 11. The nuclear metallicity evolution for four identical initial progenitor disks, set on four distinct merging trajectories, all aligned 
with final coalescence at t = 0, as indicated by the dot-dashed vertical line. Each orientation has a unique metallicity evolution, but all 
share common features. First, depressions of the nuclear metallicity are common following any pericentric passage (though not all pericenter 
passages induce nuclear metallicity depressions, as demonstrated by the "k" orientation). Second, depressions of the nuclear metallicity are 
common preceding final coalescence. 
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Figure 12. The average metallicity evolution for a set (orien- 
tations a-p) of major mergers with identical progenitor galaxies 
on different merger orientations. The dark blue region denotes the 
mean nuclear metallicity plus the one sigma variation for BB galaxy 
mergers on all 16 orientations. The light blue region shows the nu- 
clear metallicity for the same mergers when chemical enrichment 
from star formation is neglected. First pericenter passage and final 
coalescence for all orientations are denoted by the vertical dashed 
and dot-dashed lines, respectively. 



pression is identified as being on a nearly retrograde ori- 
entation. More concretely, the induced gas inflow fol- 
lowing pericenter passage is found to be small compared 
to the other systems. The idea that some orientations, 
specifically those on retrograde orbits, have suppressed 
responses duri ng pericentric passages has been stud- 
ied p reviously (jToomre fc: Toomr3ll972l : iD'Onghia et al.l 
l201Clf ). and we note that the metallicity evolution is 
affected accordingly in the cases where strong post- 
pericenter passage responses are not found. 

Subsequent depressions in the nuclear metallicity, can 
be identified as occurring during periods of strong gas in- 
flow caused by close encounters. For the "e" and "h" ori- 



entations, only one additional close encounter occurs be- 
fore proceeding to coalescence. However, for the "f" ori- 
entation, three metallicity dips can be seen. These metal- 
licity dips are natural products of this particular merger 
orientation, which leads to three distinct pericenter pas- 
sages before the galaxies ultimately coalesce. While a 
hard and fast rule governing the metallicity depression's 
dependence on merger orientation cannot be given, two 
statements can be made: most orientations undergo a 
metallicity depression following pericentric passages (un- 
less the tidal interaction is sufficiently suppressed), and 
most orientations show metallicity depressions immedi- 
ately preceding final coalescence. Both of these state- 
ments are demonstrated in figure 1111 where three of four 
(e,f, and h; not k) orientations show post-pericentric dips 
in their nuclear metallicity and three of four (e, k, and h; 
not f ) orientations show dips in their nuclear metallicities 
preceding final coalescence. 

If we identify the merging time, tm, to be the time be- 
tween pericenter passage and final coalescence, we can 
average the metallicity evolution for several orientations 
to get an orientation averaged progression of the nuclear 
metallicity during the merger. Figure [T2] shows the mean 
nuclear metallicity evolution for all 16 orientations, with 
the solid region indicating the la variations. Despite the 
variation that occurs in the nuclear evolutionary tracks 
for individual orientations, the mean result has a well- 
defined pattern with reasonably small dispersion and a 
clear interpretation. Specifically, there are two notable 
metallicity depressions caused by gas infiows, one fol- 
lowing pericentric passage and one preceding black hole 
coalescence, with a chemical enrichment induced peak 
separating in between. Also shown in Figure[T2]is the av- 
erage dynamical metallicity. Here we can see that simu- 
lations without chemical enrichment substantially under- 
predict the resulting nuclear metallicity. However, even 
when chemical enrichment is included, the mean nuclear 
metallicity is depressed throughout the interaction. 

4.3. Progenitor Masses 

We examine the effect that system mass has on the nu- 
clear metallicity evolution using a set of major mergers 
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Figure 13. The average metallicity evolution for a set (the e, f, 
h, and k orientations) of major mergers in four mass bins. 



(i.e. 1:1 mass ratios) with progenitor galaxies of vary- 
ing mass. As detailed in §3.11 the four galaxies used 
in this paper are constructed to be self-similar. Hence, 
the mass study performed here isolates the influence of 
merger mass from other parameters that may scale with 
mass (e.g. gas fraction, galaxy structure, etc.). 

We use four orientations in each mass bin (the "e", 
"f " , "h" , and "k" orientations) , and average the result- 
ing metallicity evolutionary tracks. The result, shown in 
Figure [131 is that the shape of the evolutionary track is 
nearly preserved, while the magnitude of the depression 
depends weakly on mass. We note that the preserved 
shape of the metallicity evolution supports our argument 
in the previous sections that the metallicity evolutionary 
track is a natural byproduct of the generic merger pro- 
cess. This indicates that, within the limited range of 
galaxy masses that are used here, the system mass is not 
a major consideration in determining the nuclear metal- 
licity evolution. 

4.4. Gas Fraction 

Initial gas fraction is an important parameter to con- 
sider in nuclear metallicity evolution. For higher gas frac- 
tion galaxies that, by definition, have lower initial stellar 
mass, two basic changes occur: The initial nuclear metal- 
licity of the system, as prescribed by the MZ relation, 
is lowered (decreasing the importance of the dynamical 
metallicity), and higher star formation rates occur (in- 
creasing the importance of chemical enrichment). The 
combined result, as discussed in this section, is that the 
nuclear metallicities of interacting systems at higher gas 
fractions tend to evolve to larger values, rather than de- 
creasing, in contrast to their lower gas fraction counter- 
parts. 

We first consider the nuclear metallicity evolution for 
a series of disks with varied initial gas fractions (8, 20, 
40, and 60%). We merge together these systems on the 
"e", "f", "h", and "k" orientations. If we initialize the 
metallicity properties of these disks using the local MZ 
relation, an assumption that we will relax later in this 



section, we can track the metallicity evolution and plot 
the change in the nuclear metallicity, as in Figure [141 
This shows that the depressions in the nuclear metallicity 
that occurred for the low initial gas fractions considered 
in the previous sections shift toward less negative, or 
even positive values, as the initial gas fraction is raised. 
This transition is caused by the increasing dominance of 
the enriched metallicity over the dynamical metallicity, 
primarily owing to the higher star formation rates in the 
higher gas fraction systems. 

These highest gas fraction systems, however, are not 
natural analogs to galaxies in the local Universe. There- 
fore, our choice to initialize the metallicity properties 
of these systems using the local MZ relation comes into 
question. Instead, we should consider an MZ relation 
better suited for higher redshifts. A full theoretical un- 
derstanding of the origin of the MZ relation has not yet 
been obtained, however, cosmological simulations indi- 
cate that the shape of the MZ relation is achieved at early 
times (i.e. z>6) and that galaxies move along this MZ 
relation toward higher stellar mass, while th e amplitude 



of the MZ relation slowly increases over time (iD ave et al 



in broad consistency with observations f Er b et al 



2006at iMaiolino et al.l 120081 ). Without committing our 



selves to a single choice for the high redshift MZ rela- 
tion, we can re-initialize the metallicity properties of our 
galaxies using an MZ relation that is offset by a constant 
fixed value, (5, toward lower metallicities. 

Since our simulations have no metallicity dependence, 
re-normalizing the MZ relation has no effect on the en- 
riched metallicity, but lowers the impact of the dynamical 
metallicity. The result, shown in the right panel of Fig- 
ure [141 is that as the the MZ relation is shifted to lower 
metallicities, the change in the nuclear metallicity is ex- 
pected to become increasingly positive. For all choices 
for the initial MZ relation, the higher gas fraction sys- 
tems have a median change in metallicity that is more 
positive than their low gas fraction counterparts. 

Using the local MZ relation (i.e. the black line in Fig- 
ure [14]) for systems with relatively low gas fractions, we 
find that negative changes to the nuclear metallicity are 
expected in a merger. This is consistent with the ob- 
served lowering of the MZ relation. However, as we move 
to higher gas fractions, or MZ relations better suited for 
higher redshifts, we find that increases to the nuclear 
metallicity are expected. This result is caused by shift- 
ing the dominant driver of the nuclear metallicity from 
metallicity dilution to chemical enrichment. Moreover, 
if the gas content is sufficiently high, the galactic stellar 
component cannot efficiently torque do wn the gas, re- 
sulting in moderated early stage inflows (jHopkins et al.l 
2009b). 

The same merger fundamentals apply to low and high 
gas fraction disks, independent of their initial metallici- 
ties, but the relative increase in the role of chemical en- 
richment with respect to the dynamical metallicity causes 
higher redshift or gas fraction systems to increase their 
nuclear metallicities during the merger process. 

Without being overly restrictive about our choice for 
initial gas fraction or offset from the local MZ relation, 
high-redshift, gas-rich galaxies are expected to have me- 
dian increases in their nuclear metallicity when undergo- 
ing a merger. Based on this evidence, we predict that the 
depression that occurs in the MZ relation during galaxy 



Nuclear Metallicity of Interacting Galaxies 



15 



0.20 



£ 0.15 



~ 0.10 



0.05 - 



0.00 



I 1 1 I I I I 1 1 1 I I I I 1 1 1 I I I I 1 1 I I I I 1 1 1 I I I 1 1 1 I I I I 1 1 1 I 

20% Gos Froction 

^^^^^ 40% Gos Froction 

80% Gos Froction 




-0.1 
A log(0/H) 




40 60 

Initial Gos Fraction 
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mergers and interactions in the local Universe would not 
occur at higher redshift. Instead, increases in the MZ 
relation are expected. 

5. COMPARISON WITH OBSERVATIONS 

Much work in recent years has gone into observation- 
ally classifying the metallicity of interacting o r close pair 
galaxies in the nearby Universe ([K cwIcv et al.l 120061: 
Ellison et al.ll2008l : IPeeples et al.l[2009t iSol Alonso et al I 



20101 ). In this section we compare our results directly 
with these observations. To facilitate this comparison, 
we parameterize the evolution of the nuclear metallicity 
in terms of the galactic projected separation. We define 
the projected separation as s = rcos(f) with ^ being 
a uniformly distributed random number between and 
27r. Since all of our systems are, by definition, interact- 
ing systems, this projected separation neglects contribu- 
tions from truly non-interacting systems. We present the 
mass-metallicity (MZ) and separation-metallicity (SZ) 
relations for a simulated population of interacting galax- 
ies. 

For these comparisons, we generate a population of 
progenitor galaxies that con tain the metallicity proper- 
ties of observed field galaxies (iShieldslllQQOtfBellev fc Rovl 
[1991 IZariFskv et aLl 11991 IMacArthur et al.l 120041 ). We 
use a Gaussian distribution for the nuclear metallicity se- 
le ction with a mean value consistent with the MZ relation 
of iTremonti et aLl ()2004[ ) with a constant standard devi- 
ation across mass bins of a = 0.1 dex. The metallicity 
scale length is taken from a Ga ussian distribu t ion w ith 
a mean gradient consistent with iZaritskv et al.l ()1994D of 
h;, = h/0.2 with a = 0.3h,. 

The initial metallicity properties of each galaxy in each 
snapshot are assigned based on the above criteria. Our 
ability to pick metallicity properties of the progenitor 
disks in our post-processing analysis is a consequence of 
not having a dynamical dependence on metals or allow- 
ing mixing. This lets us sample the parameter space of 
initial disk metallicity properties thoroughly, without re- 



running simulations. 

5.1. The Mass- Metallicity Relation for Interacting 
Galaxies 

The MZ relation for our simulated interacting galax- 
ies is shown in Figure [131 We select galaxies with pro- 
je cted separations less than 30 kpc - the same as was used 
in lEhisonet al.l (|2008D . We plot the MZ relation in Fig- 
ure [151 The result is a depressed MZ relation with an av- 
erage depression of 0.07 dex, which is in ag reement with 
the o bserved depression of 0.05-0.10 dex ([Ellison et al.l 
[2001) . 

Unlike luminosity, which can increase dramatically 
during merger induced starburst activity, the stellar 
mass of an interacting galaxy will change only by a 
modest factor. Therefore, the evolution of galaxies on 
the MZ relation is driven by changes to the nuclear 
metallicity, not stellar mass. In the simulation anal- 
ysis, we are able to quickly determine a system's stel- 
lar mass. The s tellar masses f or the close pair systems 
included in the lEllison et al.l (|2008[ ) sample have had 
thei r stellar masses determined via optical and IR col- 
ors ([KauflFmann et al.ll2003t ITremonti eralll2004| ) which 
have been shown to com pare well with the spectrally- 
determined stellar masses ([Drorv et al.ll200^ . Assuming 
that the observationally determined stellar masses are 
not strongly affected by ongoing starburst activity, the 
depression in both the observed and simulated MZ rela- 
tion is being caused entirely by changes to the nuclear 
metallicity - not by changes in the galactic luminosity or 
stellar mass. 

The depression in our simulated MZ relation would be 
slightly smaller if we included a set of field galaxies that 
lie on the MZ relation, but appear close in projected sep- 
aration (i.e. if we included a population of interlopers) . 
This result gives a comprehensive reproduction of the de- 
pression in the MZ relation for interacting galaxies that 
accounts for all relevant physics and the intrinsic scat- 
ter in the initial metallicity properties of the progenitor 
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Figure 15. The simulated mass-metallicity relationship for inter- 
acting galaxies is shown. The progenitor galaxy MZ relationship 
from Trcmonti ct al. ( 2004) is indicated along with the second or- 
der best fit from our simulation results. 



galaxies, while covering a reasonable portion of merger 
phase space. 

5.2. Separation vs Metallicity 

FigurefTBlshows an example of two galaxies ( "C" galax- 
ies) from one merger simulation (the "e" orientation) 
moving in separation-metallicity (SZ) space as they ap- 
proach coalescence. The evolution can be broken down 
into four distinct stages, as marked by arrows in Fig- 
ure [161 First, the galaxies approach first pericenter pas- 
sage at relatively constant nuclear metallicities. Second, 
the galaxies separate as gas inflows dilute their nuclear 
metallicity. Third, the galaxies begin falling back to- 
gether, with star formation causing a significant amount 
of chemical enrichment. And, finally, the galaxies un- 
dergo strong gas inflows that further dilute their nuclear 
metallicities as they approach coalescence. 

The details of this SZ evolution will change substan- 
tially with the merger orientation. However, we can 
gather from Figure [11] that all four stages of this evo- 
lution will be somewhat generic. Namely, galaxies will 
present constant nuclear metallicities as they approach 
each other for the first time, and strongly depressed 
nuclear metallicities as they approach final coalescence, 
with mildly depressed nuclear metallicities in between. 

Using the evolutionary track of 1161 we can develop 
some expectations for the distribution of galaxies in SZ 
space. Clearly, we expect close pair galaxies to show 
depressed nuclear metallicities. In addition, we can ex- 
pect depressed metallicities out to large galactic separa- 
tions. While the largest nuclear metallicity depressions 
are expected around close separations, widely separated 
interacting galaxies show mildly depressed nuclear metal- 
licities after interacting with their companion. The ob- 
servability of these trends can be assessed by using our 
entire merger library and converting real-separations into 
projected separations. 

The distribution of all simulated galaxies in SZ space 
is shown in Figure FTTI along with data from .Kewlev et aLl 
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Figure 16. The evolution of two galaxies from one merger sim- 
ulation are shown. The numbered arrows follow the evolution of 
the galaxies as they: 1) Approach first pericenter passage with an 
unchanging nuclear metallicity, 2) increase their separation while 
diluting the nuclear metallicity, 3) come together for the second 
time with notable contributions to the nuclear metallicity from star 
formation, and 4) lower their nuclear metallicity as they approach 
final coalescence. 



(|2006f) . The distribution of our simulated galaxies, shown 
with the colored contours in the lower half of Figure [T71 
coincide with many of the observed systems, which is 
demonstrated in the histogram in the top half of Fig- 
ure 1171 There is a trend toward lower metallicities with 
close separations. However, this trend is very mild and 
the most-extreme low-separation low-metallicity objects 
from the observed sample are not reproduced by our 
models. Despite this difference, it should be noted that 
our simulati ons agree with the m ajority of the observed 
objects from lKewlev et aLJ ()2006D . and, there is a system- 
atic but subtle drop in the simulated nuclear metallicity 
as the galaxies decrease their projected separation. 

We note that the overall distribution of simulated 
galaxies in SZ space is relatively featureless, as evidenced 
by the contours in Figure 1171 The nuclear metallicity 
depressions present at both small and large separations 
in Figure \W\ are washed out as a result of projecting 
the physical separations, the wide variety of metallic- 
ity dilution versus separations that result from different 
spin-orbit couplings, and the scatter in the initial nu- 
clear metallicities as prescribed by the mass-metallicity 
relation. In addition, while late stage mergers (e.g., stage 
4 in Figure [TC)) present the largest depressions in nuclear 
metallicity, they only last for a short period of time com- 
pared to the overall merger process. Since we place equal 
weight on every timestep for every simulation when de- 
termining the distribution of simulated galaxies in Fig- 
ure [T71 the highly-depressed late-stage mergers will nat- 
urally have a less dramatic impact. 

Several factors can influence the distribution of our 
model galaxies in SZ space. In particular, the 
(in)completeness of our merger library and the metal en- 
richment scheme will influence the resulting SZ relation. 
We discuss both of these factors in the following subsec- 
tions, so that we may better understand the disagreement 
between our models and the observations at low separa- 
tions and low metallicities. However, we also note that 
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Figure 17. The relative frequency of galaxies in metallicity space 
(top) and the sim ulated separat ion-metallicity relationship along 
with the data from lKewlev et alj (2006) (bottom) arc shown. The 
majority of observed points from Kcwlcv et al. (2006) are coinci- 
dent with our simulated galaxies, with the exception of the low- 
metallicity, low-separation tail. Our simulated galaxies show a de- 
pression in the nuclear metallicity for close pair galaxies, however 
this depression is less dramatic than is required to explain the low- 
metallicity, low-separation tail observed iu iKewlev et al.i t2006i '). 



additional observations would help clarify how common 
these low-separation, low-metallicity objects are. 

5.2.1. Merger Library Completeness 

Although the merger library used in this work was 
fairly comprehensive, there is a significant portion of 
merger parameter space that has not been explored. Per- 
haps this could explain why we do not reproduce the low- 
separation, low-metallicity objects. We can describe our 
exploration of merger parameter space with two compo- 
nents: Variations of initial isolated galaxies models and 
merger configuration parameters. 

Additional merger configurations (i.e. variations to 
the orbital orientation and angular momentum) are un- 
likely to yield surprising results. Out of all of our 
merger simulations, a maximum depression of about 0.3 
dex was found, with many depressions being must less 
dr astic. To achieve the lowest-metallicity objects found 
in IKewlev et al.l ()2006l ). depressions of ~ 0.6-0.7 dex 
would be required. It seems unlikely that the isolated 
galaxy models used in this study will reproduce very low 
metallicity objects, even with a more extensive explo- 
ration of merger orientations. 

On the other hand, galaxy models that are signifi- 
cantly different from the ones used in this study, such 
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Figure 18. The simulated separation-metallicity relationship is 
shown using a delayed chemical enrichment scheme. Newly formed 
metals are allowed to contribute to the HII region metallicity mea- 
surements only after the SPH particle in which they reside under- 
goes an integrated star formation rate equal to its own mass. This 
effectively delays metal enrichment by one generation. 



as low mass galaxies, could fill out a new region in SZ 
space. For example, galaxies with initial stellar masses 
of Ms ^ lO^Af0 will have initial nuclear metallicities of 
12-|-log(0/i?) ~ 8.6 and could have their nuclear metal- 
hcities depressed to ~ 8.3 or 8.4 during the merger. More 
generall y, it seems lik e ly tha t the lowest metallicity ob- 
jects of IKewlev et al.l (|2006( ) are likely merging systems 
that had low metallicities prior to the onset of the merger 
which may be a reasonable possibility given the low B- 
band luminosities associated with these objects. These 
systems may be reproduced by our simulations when we 
include a more extensive set of initial galaxy models. 

5.2.2. Approximating Delayed Mixing 

Although oxygen abundances are well-approximated 
by an instantaneous enrichment model, this does not 
necessarily imply that metals produced in SNe instan- 
taneously migrate to HII regions. Specifically, we can 
relax our assumption that the metallicity of an SPH par- 
ticle should be instantaneously identical to the metallic- 
ity of an HII region. Our analysis calculated the nuclear 
metallicity from the current metallicities of gas particles. 
While the metals produced in core collapse supernovae 
are quickly returned to the ISM (~ 10^ years) and su- 
pernovae remnants cool and fade into the ambient ISM 
over similar timescales (~ 10^ years), there could then 
be some additional delay before this gas is sufficiently 
mixed with the surrounding medium and ends up in HII 
regions. In particular, we expect the metals produced by 
SNe to affect the next generation of HII regions. As a 
first approximation, this does not require a change to our 
simulations, but rather a change to our post-processing 
analysis. 

We can define an effective HII region metallicity that 
is given by the metallicity that the SPH particle had 
one generation of star formation ago. If the number of 
generations of star formation is given by 



1 r dhU , 



(17) 
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Figure 19. The simulated separation-metallicity relationship is 
shown using a delayed chemical enrichment scheme. Newly formed 
metals are allowed to contribute to the HII region metallicity mea- 
surements after a specified delay time, t^j, in Gyrs. 



where M is the particle mass, then the effective HII re- 
gion metal content is given by 

f Mz,^n^t + yM (TV* - 1) if A^* > 1 

M'zit) = { (18) 

[ Mz^^n^t if /V, < 1 

which reflects the particle's metallicity one generation of 
star formation in the past. 

This delayed chemical enrichment scheme gives an SZ 
relation, shown in Figure [T8l which has a low metallicity 
tail that is qualitatively similar to the observed galax- 
ies. This does not imply that our previous calculation 
of the nuclear metallicity using the non-delayed chem- 
ical enrichment values were incorrect, but rather that 
the assumption of instantaneous transition of metals into 
HII regions may be overly restrictive. This is a plausi- 
ble explanation, provided that the associated delay-times 
are reasonable estimates for the timescale over which SN 
ejecta would mix and diffuse into the ISM. 

If we parameterize an effective HII region metallicity 
in terms of a delay-time, rather than an integrated star 
formation rate, we can get an estimate for the required 
delay time. Rather than neglecting the enrichment from 
the last generation of star formation, we neglect the en- 
richment that occurred over the past td years. The effec- 
tive metal content of a particle is then given by 



M'zit) = 



Af 



ZAjiit 



+ yIo 



t-td dM. 
dt 



M 



Z.init 



'i^dt' if i > td 



if i < td 



(19) 



where td is the assumed time delay. 

This time delayed chemical enrichment scheme gives an 
SZ relation, shown in Figure [TOl where larger delay times 
lead to systematic decreases in the inferred metallicity of 
close pairs. This is similar to those seen in Figure [181 but 
with an explicit time dependence. The delay times re- 
quired for our simulations to match observations is of the 
order of a few 10^ years, which is comparable to estimates 
for the metal mixing timescale of newly produced met- 
als into the ISM (,Oev,,200a : ,Scalo fc Elmegreen.,2004, ). 



We conclude that delaying the transfer of metals from 
SN remnants to HII regions is important to consider; 
however, it is unclear whether these long mixing times 
are appropriate for the extreme environments in merg- 
ing galaxy nuclei where external forces, gas inflows, high 
gas densities, and large star formation rates are present. 

6. DISCUSSION 
6.1. Predictions for Galaxies with High Gas Content 

The results of § 14.41 demonstrate that the merger- 
induced metallicity evolution is heavily dependent upon 
gas content, specifically, higher gas fraction disks pro- 
duce less nuclear metallicity depression than lower gas 
fraction disks (see, e.g.. Fig. [Tl)) . In the extreme, high 
gas content mergers may even yield significant enhance- 
ments of central metallicity, a result that emerges both 
because of the increased impact of enrichment from star 
formation as well as because nuclear gas i nflows are sup- 
presse d when the gas fraction is very high (jHopkins et al.l 
l2009bD . 

As a consequence of these trends, we expect that the 
nuclear metallicities of gas-rich merging galaxies will be 
similar to, or greater than, an appropriately matched 
quiescent sample. This is in contrast to the depression of 
nuclear metallicity observed in low gas-content mergers 
which was the focus of this paper. Such considerations 
may be applicable at high redshift where current obser- 
vations suggest that galaxies contain significantly higher 
gas fractions (Tacconi et al. 2010; Erb et al. 2006b) and 
suppressed metallicity at fi xed galaxy mass (JErb et al.l 
l2006at [Maiolino et al.ll2008l) o r for less-massive gas -rich 
galaxies in the local universe (jCatinella et al.|[2010[ ). 

The predicted metallicity enhancements in gas-rich 
merging systems was also noted in a recent paper by 
iPerez et al.l (|2011| ) using similar methods. While both 
works observed similar metallicity enhancements, the 
timing was not identical, and purported cause o f these 
effects were very different. The IPerez et al.l (|2011f) model 
shows rapid metallicity enhancement during early stages 
of the merger primarily as a result of rapid star forma- 
tion owing to disk instabilities in the gas-rich system 
(see, e.g., their Figs. 8 and 9). Such effects occur prior 
to, or during the early stages of the merger and are there- 
fore decoupled from the final coalescence stage. In our 
model, the gas-rich disks are designed to be stable a 
priori, and therefore the metallicity enhancements are 
driven by intense star formation triggered by the merger. 
While the end result is the same, namely metallicity en- 
hancements, the relative timing and strength of the vig- 
orous star formation, the metallicity dilution owing to 
infiows, and the metallicity enhancement owing to star 
formation are fundamentally different, and therefore po- 
tentially discernible via observations. 

We caution, however, that the galaxy merger simula- 
tions carried out here may require some generalization 
before they are directly applicable at high redshifts. Ef- 
fects such as the higher UV background, higher merger 
rate, and direct fueling of low metallicity gas via the 
surrounding intergalactic medium at high redshifts may 
be of significant influence in the observed metallicity of 
galaxies. And numerical uncertainties, such as the proper 
treatment of feedback from star formation and accreting 
black holes will invariable need to be explored. 
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6.2. Future Considerations: An Improved Chemical 
Enrichment Model 

By comparing our models of interacting galaxies to 
the observed mass- metallicity (MZ) and separation- 
metallicity (SZ) relations we have refined our under- 
standing of the factors at play during the metallicity evo- 
lution of galaxies. It is crucial to point out (yet rarely 
done so), however, that numerical results such as ours are 
based upon the ability to robustly track a number of com- 
plicated and interrelated astrophysical processes. While 
we consider our sub-resolution models to be well moti- 
vated and well tested, we admit that there are straight- 
forward improvements to these models which will allow 
for more stringent tests of our physical picture. We now 
take a moment to outline these future improvements to 
our models. 

The chemical enrichment model employed here tracks a 
global metallicity through a single scalar variable that is 
continuously updated assuming that the entire solar yield 
(assumed to be a mass fraction of 0.02) is instantaneously 
recycled within a single star-forming fluid element (i.e., 
an SPH particle). We admittedly note the simplicity of 
this model, but consider the ease with which it is im- 
plemented, the straightforward analysis it allows, and its 
inherent self-consistency as distinct advantages of this 
approach. The instantaneous recycling approximation 
is also justified because we compare to oxygen derived 
metallicities, and oxygen is thought to be recycled on a 
short (~ 10^ years) timescale through core-collapse su- 
pernovae. Moreover, in a practical sense, similar models 
have been employed in nearly every SPH code over the 
last two decades which makes our results easy to gauge 
within the existing theoretical framework. 

There are, however, avenues in which distinct progress 
can be made to our chemical enrichment model. Ef- 
forts toward this end have been pioneered by a num- 
ber of authors over the last ^5-10 years (see , e.g., 
K awata 2001; Tornatorc et al. 200 7; Scannap i eco et al.l 
20051 : lOppenheimer fc Pavel 120081: iKim et al.l r2009l to 
name a few) by employing two common features: 
metallicity-dependent cooling, and more complex enrich- 
ment schemes that explicitly track various species and 
their time-dependent release. Our gas recycling scheme 
is particularly well-suited for including mass return by a 
single stellar population. As individual star particles are 
converted back to interstellar gas, they can straightfor- 
wardly carry with them the appropriate enrichment from 
Type I and II supernovae as well as mass-loss from mas- 
sive young stars and asymptotic giant branch stars. Such 
an approach is necessary to study trends in abundance 
ratios, such as those observed in the centers of elliptical 
galaxies. 

Additional accuracy may also be gained by includ- 
ing sub-resolution models for th e turbul ent diffusion o f 
metals (see, e.g.. Wad slev et al.| l2008: Sh en et al.|[20Tol ). 
Coupled to an advanced enrichment algorithm, a physi- 
cally motivated procedure for diffusing metals will allow 
enrichment to occur locally but still have global influ- 
ences if conditions warrant. Such methods have been 
shown to be important f or tracking metal licity in cosmo- 
logical simulations fe.g.. lShen et al.l l2010') and we intend 
to assess their influence in isolated and merging galactic 
systems. 



Finally, the comparisons between the models and the 
observations can be significantly improved by actually 
tracking the detailed h ne emission using 3D radiative 
transfer (jJonsson et al.l [2009) . As efforts to produce data 
cubes with full spatial and spectral information are de- 
veloped (e.g. Rich et al. 2010) . full metallicity maps can 
be produced and directly compared to our simulations. 
Radiative transfer in our simulations will shed light on 
any additional information that may be encoded in the 
line profiles and allow for a more even-handed compari- 
son. 

7. CONCLUSION 

Using numerical simulations, we have investigated the 
impact that mergers have on galactic nuclear metallic- 
ity. Our models include the capability of describing star 
formation, chemical enrichment, gas recycling, and star- 
formation driven winds. The analysis performed here re- 
lies upon the ability to accurately track the spatial and 
temporal return of metals to the ISM enabled by our 
stochastic gas-recycling algorithm within GADGET (see 
§[22|. 

One of the primary results of this work is to rein- 
force the notion that galaxy mergers, and their attend- 
ing gravitational tidal forces, generate significant infiows 
of gas shown exphc i tly in Figures [TT1 and [T^ fsee also 
Rupke et all l2010al: iMontuori et"aLn2010nPerez et al.l 
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nuclear region and reinforce the generic connection be- 
tween close galaxy interactions and nuclear metallicity 
dilution - no matter how many close passages occur. 

Merger-induced gaseous inflows are, however, not the 
only factor influencing the nuclear metallicity, and this 
paper quantifies the competing effects of enrichment from 
star formation (§ I4.1.ip . gas consumption (S I4.1.2|) . and 
galactic winds (§ I4.1.3|) . Our models provide a physical 
template to understand nuclear metallicity evolution and 
lead us to conclude that merger-induced metallicity dilu- 
tion and chemical enrichment from star formation drive 
the nuclear metallicity, while gas consumption and galac- 
tic winds play a secondary role to modulate the efficiency 
of the primary processes. 

We have made a distinct effort in this work to directly 
compare our models to observations, allowing us to val- 
idate the physical evolutionary model that we have pre- 
sented and associate observable trends with the various 
stages of the merger evolution. We find that the merger- 
induce depressions of nuclear metallicity are typically ~ 
0.07 dex, an offset similar to that observed to the mass- 
metallicity relation for close pairs in SDSS ([Ellison et al.l 
l2008f ). The models also demonstrate that central de- 
pressions of metallicity are time-dependent owing to the 
specific merger orbit, the galaxies' spin-orbit coupling, 
and varying structural properties of the galaxies. When 
co mpared to the obser ved separation-metallicity relation 
of lKewlev et al.l (|2006[ ). the models generally show a sim- 
ilar trend to produce lower central metallicities at smaller 
separations. The precise role of interactions is difficult 
to discern, however, because of the significant scatter in 
the relations (see [17] and the discussion in §5.2p . Future 
comparisons in which galaxies are morphologically clas- 
sified according to merger stage, or separated according 
to IR luminosity or star-formation rate may be better 
suited to isolated the true effects of the interaction. 
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This work also demonstrates that the progenitor gas 
content can have a profound influence on the resulting 
merger-induced metallicity dilution. Mergers between 
gas-rich progenitors can yield systematic metallicity en- 
hancement, as opposed to metallicity dilution (see i; 14.41 
and specifically Fig.[T4|). In our picture, central metallic- 
ity enrichment comes from vigorous merger-induced star 
formation that can compensate for, and eventually over- 
come the merger-induced inflow of metal-poor gas. This 
vie w i s slightly differ ent than put forth in recent work 
by [Perez ct al. ()2011| ) , which argues that disk instabili- 
ties and clump formation drive rapid star formation and 
subsequent metal enrichment. But once these processes 
occur, which is typically early in the merger process, the 
story becomes similar to when the galaxies have low gas 
content. By studying the detailed correlations between 
merger stage, pair separation, central metallicity, and 
metallicity gradients, in both observations and the mod- 
els, future work will allow for a better understanding of 
the (merger and enrichment) processes at work. 
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